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ABSTRACT 



Context. The recently launched satellite, Fermi Gamma-ray Space Telescope, is expected to find out if cosmic-ray (CR) 
protons are generated from supernova remnants (SNRs), especially RX J1713.7— 3946, by observing the GeV-to-TeV 
7-rays. The GeV emission is thought to be bright if the TeV emission is hadronic, i.e., of proton origin, while dim if 
leptonic. 

Aims. We reexamine the above view using a simple theoretical model of nonlinear acceleration of particles to calculate 

the gamma-ray spectrum of Galactic young SNRs. 

Methods. 

Results. If the nonlinear effects of CR acceleration are considered, it may be impossible to distinguish the evidence of 
proton acceleration from leptonic in the 7-ray spectrum of Galactic young SNRs like RX J1713.7— 3946. On the other 
hand, future km -class neutrino observations will likely find a clear evidence of the proton acceleration there. 

Key words. Acceleration of particles - ISM: supernova remnants - Gamma rays: theory 



1. Introduction 

Recently, Fermi Gamma-ray Space Telescope (FermiJ^, ob- 
serving GeV 7-ray photons, has been launched. The GeV 
7-ray observations with Fermi are expected to identify the 
accelerators of Galactic cosmic-ray (CR) protons whose 
energy extends up to the "knee" energy (« 10 15 5 eV). 
At present, the most probable candidate for the CR ac- 
celerator is a young supernova remnant (SNR). Since the 
detections of synchrotron X-ray s in some SNRs sho w ev- 
idence of electron acceleration (|Kovama et alJ ll995). the 
current unsolved issue is whether the SNRs produce high- 
energy protons or not. TeV 7-ray observations are impor- 
tant to address this problem. So far, TeV 7-rays have been 
detected from several young SNRs (lEnomoto et alj [20021; 
lAharonian et alj 12004 120051 ; iKatagiri et alj I2005D . They 
arise from either leptonic (CMB photons up-scattered by 
high energy electrons) or hadronic (7r°-decay photons gen- 
erated via accelerated protons) processes, and it is generally 
difficult to separate these processes using only the TeV en- 
ergy band; the study of wide-band, GeV-to-TeV spectra is 
necessary. 

RX J1713.7-3946 (hereafter RXJ1713) is a representa- 
tive SNR from which bright TeV 7-rays have been detected. 
The H.E.S.S. experiment measured the TeV spectrum and 
claime d that its shape was better ex plained by the hadronic 
model (jAharonian et al . 2006, 2007). So far, compared with 
other young SNRs, the TeV 7-ray spectrum of RXJ1713 is 
the most precisely measured and the energy coverage is 



wide, from 0.3 to 100 TeV, so that we can obtain the best 
constraints on theoretical models. 

Recently, time var iation of synchrotron X -rays was dis- 
covered in RXJ1713 (jUchivama et al.ll2007h . If the varia- 
tion timescale is determined from the synchrotron cool- 
ing of X-ray emitting electrons, the magnetic field is es- 
timated to be B ~ mG. If so, the leptonic. one-zone emis- 
sion model (e.g. JAharonian fc Atovanll999f ) cannot explain 
the TeV-to-X-ray flux ratio, supporting the hadronic ori- 
gin of TeV 7-rays. It should be noted that the amplified 
magnetic field is theoretically ex pected (e.g. jLucek fc Bell 
I2000t iGiacalone fe Jokipiil [2007h . In this case, according 
to the standard diffusive shock acceleration theory, the 
maximum energy of accele rated protons is estimated as 
(jAharonian fc Atovanlll999h 



-Emax.p — 8 X 10 3 



Vg 



.4000 km s- 



TeV 



(1) 



http:/ /fermi. gsfc.nasa.gov/ 



which can be comparable to the knee energy. Here, B m Q, 
V s , £3, and rj g are the magnetic field strength in units of mG, 
the shock velocity, the age of the SNR in units of 10 3 yr, 
and the gyrofactor, respectively. 

However, at present, there are several issues to be ad- 
dressed, as the above picture on RXJ17I3 is not yet proved. 
First, if B ~ mG and TeV emission is hadronic, then in 
order to explain the measured flux of radio synchrotron 
emitted by primary electrons, the electron-to-proton ra- 
tio at the SNR should be anomalo usly small, K ep ~ I0~ 6 
(jUchivama et alj|2003t lButtll2008f) . which is far below the 
observed value at the earth and estim ated values in the 
nearby galaxy (jKatz fc Wax man 2008). This might be re- 
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solved if the electrons are accelerated in the later stages of 
SNR evolutio n, when the value of K ep is different from the 
present value (jTanaka et al.ll2008l ), although further discus- 
sions are necessary. Second, the hadronic scenario may be 
inconsistent with the molecular cloud (MC) observations 
(|Fukui et al.ll2003D . RXJ1713 is surrounded by MCs, which 
might suggest collision with them and high target number 
density. If the TeV 7-rays a re hadronic, such a region should 
be brighter than observed {Plaga 20(51). 

In addition, if the measured width of the synchrotron 
X-ray filaments at the shock front of SNRs is deter- 
mined by the synchrotron cooling effect (lUchiyama et al.l 
120031: IVink fc Laniinel 120031 : iBamba et al.1 120031 l2005aU bT). 
the mag netic field is indep endently estimated as B w 
0.1 mG (jParizot et al.l [200a ), which is a n order of mag- 
nitude smaller than that estimated by lUchiyama et al.l 
(|2007h . Also, the cutoff energy of TeV 7-ray spectrum is 
low, so that in the one-z one hadronic scenario E max v is 
estimated as 30-100 TeV (jVillante fc Vissani|[2007l) . which 
is approximately two orders of magnitude lower than the 
knee energy. If £ lnaXiP < 100 TeV and B w 1 mG, then 
Eq. (1) tells us r\ g J; 80, implying far from the "Bohm 
limit" (rig ~ 1) which is inferred from the X -ray observation 
(|Parizot et al. 2006; Y amazaki et al.ll2004 ) or expected the - 



oretically (jLucek fc BellH2000l : iGiacalone fc Jokipiil I2007T) . 
This statement is recast if we involve recent results of X- 
ray observations. The precise X-ray spectrum of RXJ1713 
is revealed, which gives v s = 3.3 x lO 8 ^ 2 cm s _1 
(jTanaka et al.l [20081 ) . Then, Eq. (1) can be rewritten as 
-Bmax.p = 5 x 10 3 -B m Gi3 TeV. Hence, in order to obtain 
-Bmax.p < 100 TeV, we need B 20 /J.G in the context of 
the hadronic scenario of TeV 7-rays. One might think that 
the volume filling factor of the region with B « 1 mG is 
small and that the average field strength is smaller, e.g., 
B « 0.1 mG. However, even in this case, -E m ax,p is more 
than 100 TeV, which contradicts the observed 7-ray spec- 
trum beyond 10 TeV. 

In these circumstances, Fermi will give us important 
information on the 7-ray emission mechanism. So far, the 
GeV emission has been thought to be bright if the TeV 
emission is hadronic, while dim if leptonic. However, this 
argument is not so straightforward if the nonlinear model of 
CR acceleration is considered. In the next section, we cal- 
culate the photon spectrum using a simple semi-analytic 
model taking into account nonlinear effects. Indeed, we 
show that in a certain case, the hadronic emission spec- 
trum in the GeV-to-TeV band is similar to the leptonic 
one. 

2. Hadronic gamma-rays in the efficient 
acceleration case 

If a large amount of protons are accelerated, their momen- 
tum flux is large, so that the back-reaction of them is sig- 
nificant and the background shock structure is modified 
(iDrurvi [19831 iBlandford fc Eichler! [19871 : iMalkov fc Dr"u7vi 
|2001| ). Compared with the test-particle (inefficient acceler- 
ation) case in which the back-reaction effects are neglected, 
the background plasma is more compressed at the shock due 
to the additional CR pressure, which leads to a harder CR 
spectrum. Hence the hadronic emission becomes hardeiQ- 



At present, there is no reliable theory to determine the 
acceleration efficiency, and it is not clear whether this non- 
linear model is correct or not. Thus, the observations to 
determine the acceleration efficiency and the CR spectrum 
at the acceleration site are important. It is widely expected 
that RXJ1713 with precise studies in the 7-ray and X-ray 
bands is one of the best laboratories to investigate theories 
of nonlinear acceleration. 

There are several models of nonlinear CR acceleration 
(iBerezhko et al.lll994tlEllison et al.lll996l:lKang et al.ll200ll: 
i Blasil 120021: iBlasi et al.ll2005HMalkovlll997t lAmato fc Blasil 
120051) . Here, we adopt the one-dime nsional, semi-analytic 
model (|Biasill2002t IBlasi et al.ll2005f ). In the following, we 
briefly summarize the formalism. The accelerated CR pro- 
tons are described by the distribution function, f(x,p), 
where x is the spatial coordinate and p is the momentum 
of the accelerated proton. We derive stationary solutions to 
the set of an equation for f{x,p) that describes the diffusive 
transport equation of accelerated protons, and equations for 
the background thermal plasma that is treated as a fluid. 
The velocity, density, and thermodynamic properties of the 
fluid can be determined by the mass and momentum con- 
servation equations, with the inclusion of the CR pressure 
calculated as 



4^ /"Pmax 

Pcr(x) = — p 3 v(p)f(x,p)dp 

6 Jpinj 



(2) 



The injection of accelerated particles is assumed to occur at 
the shock front (a: = 0), and mono-energetic injection with 
the injection momentum p; n j is adopted. This is because at 
present we do not understand the injection process from 
first principles, hence we take a simple injection recipe. We 
assume pi n j = £pth, where pth = (2m p fcBT2) 1//2 is the mo- 
mentum of particles in the thermal peak of the Maxwellian 
distribution of the background plasma in the downstream 
region, having temperature T^. Furthermore, we assume the 
continuity of the distribution function of accelerated par- 
ticles and the background thermal plasma at p = p- m j and 

x = 0, namely /(0,pi n j) = (n 2 /7r 3 / 2 p 3 h )e _(i>ini / pth ) 2 , where 
ri2 is the downstream number density of the background 
fluid. In this formalism, input parameters are the maxi- 
mum momentum p ma x, the upstream Mach number Mq, 
upstream fluid velocity uq, and the injection parameter £. 
Given these parameters, the distribution function of accel- 
erated protons at the shock front fo(p) is calculated, and 
various physical quantities are obtained such as the injec- 
tion rate 77, the total compression ratio Rtot, the fraction of 
the CR pressure at the shock £ c (0) = -Pcr(0)//Oo u o- 

Note that recent nonlinear models of CR acceleration 
have been dev eloped taking into accou nt the magnetic field 
amplification (Vladimi rov et al]|2006|) and its influence on 
turbu lent heating (Vladimirov et al. 2008^, fluid c ompres- 
sion (|Terasawa et al.l 120071: ICaprioli et al 



[2008al[bh . and 

Alfvenic drift (jZirakashvili fc Ptuskinll2008t ). which are ne- 



glected in the model considered in this paper. These effects 
lead to less spectral hardening of accelerated particles and 
smaller compression ratios, and migh t be important in or - 
der to calculate the 7-ray spectrum ([Morlino et al.l l2008). 
However, at present, it is not certain whether or not the 



2 Another kind of formation of a hard 7-ray spectrum from ac- 
celerated protons is the SNR-MC interaction system with appro- 



priate separation (e.g.. lAharonian fc At ovan 1996; Gabic fet al.l 
l2007f ). The slower propagation of the low-energy protons toward 
the cloud makes the 7-ray spectrum hard. 
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Fig. 1. ^F^-spectra in 7-ray band, predicted by hadronic 
7r°-decay model in the efficient acceleration case (solid line) 
and leptonic IC model in the inefficient case (dashed line). 
The dot-dashed curve shows the f year, 5 a sensitivity for 
the Fermi L AT taking into acco unt the Galactic diffuse 
background (iHigashi et al.ll2008f). The observed spectrum 
in the TeV band is shown (jAharonian et al.l l2007) . 

magnetic field is strongly amplified in the acceleration re- 
gion; although streaming instabilities between accelerated 
protons and background plasma may occur, the nonlin- 
ear evolution of the instability and the saturation level 
are highly uncertain. Although the magnetic field ampli- 
fication is potentially coupled to the high injection rate of 
protons, they should, in principle, be treated separately. In 
this sense, our model is the extreme limit of the nonlinear 
acceleration theory, which predicts the hardest spectrum of 
accelerated particles. 

In this paper, we adopt Mq = 100, p max = 1 X 10 m p c, 
uq = 5 x 10 8 cm s _1 , and £ = 3.6. Then, we obtain 
p ini = 3.74 x 10~ 3 m p c, 77 = 2.03 x 10~ 4 , R to t = 36.2, 
and £ c (0) = 0.902. While the total number of CR pro- 
tons is much smaller than that of the bac kground plasma 
(77 <C 1), the CR pressure is dominant (jAmato &: Blasil 
2005). We find that the CR energy spectrum is asymptoti- 
cally N p (p) oc p 2 fo(p) ocp" 1,5 , which is harder than in the 
case of inefficient acceleration. This result on the asymp- 
totic form has been analytically derived, which does not 
depend on the shock parame t ers su ch as Mq and uq in the 
large-Mo limit (|Malkovlll99"7lll999D . 

Using the derived distribution function of CR protons, 
we calculate the 7-ray spectrum produced by 7r°-decay 
proce ss. We used the PYT HIA Monte-Carlo event gener- 
ator (Si ostrand et al.l l2006), which fits existing experimen- 
tal data well, to calculate the pp scattering processes and 
detailed distributions of the daughter particles such as ir° 
and 7r ± . We have also obtained the distribution functions 
of emitted photons and neutrinos which are produced by 
subse quent decays of those mesons and muons in the same 
code (|Yamazaki et al.| [2006'). The result is shown in Fig. 1. 

So far, we have considered the hadronic 7-ray spectrum 
in the context of the efficient acceleration scenario. For 
comparison, we show, in Fig. 1, the spectrum of leptonic 
inverse-Compton (IC) radiation via accelerated electrons 



in the case of inefficient acceleration, where back reaction 
effects of accelerated protons are neglected. Throughout the 
paper, we take into account the Klein-Nishina effect in cal- 
culating the IC spectrum. The assumed form of the elec- 
tron distribution is N e (E e ) oc E~ Se exp(— E e /E max ), and 
we adopt s e — 2.0 and E maK = 28 TeV. This case can be 
realized if the magnetic field is weak enough for the syn- 
chrotron cooling effect to be insufficient, whose condition 
is written as i agc < i S ynch = 67tto 2 c 3 /arE max B 2 , where 
*agc = £3 x 10 3 yr is the age of the SNR. Solving this equa- 
tion with E max = 28 TeV, we derive B < 26t^ 1/2 /aG. If 

the magnetic field is strong (B ^> 26i 3 /jG), the spec- 
tral deformation occurs, which will be discussed in § 0J 
One can find from Fig. 1 that the 7r°-decay 7-ray emission 
in the efficient acceleration model coincides with the lep- 
tonic IC model in the inefficient case. The reason is simple. 
Let Si be the index of the energy spectrum of accelerated 
particles i (i = p or e), so that Ni(Ei) oc E^ s * . Then, 
the radiation spectrum of 7r°-decay 7-rays is in the form 
vF v oc v 2 ~ Sp , while the spectrum of IC radiation is given by 
vF v oc jy _ ( s '=~ 3 )/ 2 . Hence, hadronic emission with s p f=a 1.5 
and IC emission with s e « 2.0 give the same spectral slope. 
This is summarized in Table [TJ 

Below several hundreds of MeV, hadronic 7-ray emis- 
sion is dimmer than leptonic IC emission because 7r° cre- 
ation reaction does not occur for low-energy (< 70 MeV 
in the center-of-mass frame) protons. Unfortunately, Fermi 
sensitivity is not high enough to recognize this decline be- 
low - GeV. 

One can find that both models slightly deviate from the 
observed spectrum in the sub- TeV energy range. The signif- 
icance of this is sometimes strengt hened, because the lep - 
tonic one-zone IC model is unlikely (jAharonian et all2 007). 
However, as will be seen in the next section, it is not serious 
if the two-zone models are considered. 

3. Two-zone models 

Here, we consider simple two- zone models to better explai n 
the observed TeV spectrum (jAharonian &: Atovadll99ll . 
RXJ1713 is interacting with MCs, so that the environment 
around the shock producing high-energy particles may be 
inhomogeneous. In this case, the one-zone approximation is 
too simplified, which motivates us to investigate the two- 
zone model as the next-order approximation. 

3.1. Hadronic two-zone model 

In this model, two independent regions, j (j = 1,2), are 
considered. For each component, we independently calcu- 
late the proton spectrum again using the semi-analytic 
model of nonlinear CR acceleration considered in the pre- 
vious section. The region j has parameters Mq , p max ^, 
Uq , and £U) . Then, we derive the hadronic 7-ray spectrum 
produced by 7r°-decay process. The total emission spectrum 
from the SNR is simply given by the sum of the emissions 
from two regions. 

Figure 2 shows the result where we adopt p max ^ 1 ' > = 
2 x 10 4 to p c and p max ^ 2 ^ — 2 x 10 5 m p c. The rest of the 
parameters are the same as those of the previous section: 
M^ l) = M (2) = 100, 4 1} = 4 2) = 5 x 10 8 cm s _1 , and 
£W = £;( 2 ) = 3.6. The normalization of the 7r°-decay emis- 
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Table 1. Spectral index of vF v spectrum of gamma-rays (yF v oc v a ) for various cases. Hadronic emission model in the 
case of efficient acceleration (la) predicts similar 7-ray spectral slope with the leptonic, inefficient acceleration model 
with weak magnetic field (Illb) . On the other hand, hadronic inefficient acceleration models (Ila and Ilia) predict similar 
7-ray spectral slope with the leptonic, moderate magnetic field model (lib). 





(a) 7r u model 
a ~ 2 — s p 


(b) IC model 

a = -(s e - 3)/2 


I Efficient acc. (strong B-field) 


(la) s p w 1.5, a w 0.5 


(lb) s e « 2.5, a w 0.25 


II Inefficient acc. (moderate B-field) 


(Ila) s p « 2, a w 


(lib) s e K3,tt!s0 


III Inefficient acc. (weak £>-ficld) 


(Ilia) s p « 2, a » 


(Illb) s e ^2,afs 0.5 
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Fig. 2. j/F^-spectra in the 7-ray band, predicted by the 
two-zone hadronic model (solid line). The dashed and dot- 
ted lines represent fluxes from region 1 and 2, respectively. 
Others are the same as in Fig. 1. 
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Fig. 3. ^-Fy-spectra in the 7-ray band, predicted by a two- 
zone leptonic IC model (solid line). The dashed and dot- 
ted lines represent fluxes from region 1 and 2, respectively. 
Others are the same as in Fig. 1. 



sion is proportional to the product of the amount of the ac- 
celerated protons, which is represented by N p (p = mc), and 

the target number density, n t . Here, we adjust rij. Np (p = 

(2) (2) 

mc)/n t Np (p = mc) = 0.56 in order to explain the ob- 
served 7-ray spectrum. Then, one can see that the fit be- 
comes better compared with the one-zone hadronic model. 

3.2. Leptonic two-zone model 

The observed c orrelation between TeV 7-ray and syn- 
chrotron X-rays (Ah aroman et al.ll2006l) may suggest that 
they have the same origin. Then, since synchrotron X-rays 
arise from accelerated e lectrons, one may expect that the 
leptonic model is likely ( | Lazendic et alj 2004; Po rter et alj 
I2006t lOeasawara et al.ll2007fl . 

Similar to the hadronic two-zone model, two indepen- 
dent regions are considered. The region j (J — 1,2) has a 
magnetic field B^> and the electron spectrum N^(E e ) — 

A^Ee S ° exp(—E e /Emlx), where is the normaliza- 
tion constant. Here, the electron spectra are given by a 
single power-law form, because the inefficient acceleration 
is adopted. We consider synchrotron emission and IC emis- 
sion in which the target photon is the CMB. The total 
emission spectrum from the SNR is given by the sum of 
the emissions from two regions. 



Figures 3 and 4 show the result where we adopt B^ — 
2.1 fiG, BW = 10 fj,G, si 1} = si 2) = 2.0, E^L = 10 TeV, 
E^L = 40 TeV, and A^/A^ = 2.74. The observed spec- 
trum, including radio and X-ray bands, can be explained 
by this model. Note that in the leptonic model, the mag- 
netic field strength m ust be much less than the observation- 
ally inferred values (jUchivama et alj l2007t iParizot et alj 
I2006D in order to fit the radio and X-ray synchrotron spec- 
trum — if the magnetic field were larger than 10 /iG, the 
predicted sync hrotron radiation would be much brighter 
than observed (|Aharonian fc Atovanl 119991 ). Hence, other 
explanations for the observations of rapid time variability 
and thin width of synchrotron filaments may be nec essary 
(|Pohl et al.ll2005t lBuTl[20Cjl iKatz fc Waxmanll2008f) . 

Comparing Fig. 2 with Fig. 3, we find that the predicted 
spectrum of the hadronic two-zone model in the efficient 
acceleration case is similar to that of the leptonic two-zone 
model in the case of inefficient acceleration. This conclusion 
is the same as that for one-zone models in § 2. 

4. Conclusion 

Some models predict relatively bright GeV 7-rays compared 
with those considered above. If the CR back-reaction effect 
on the particle spectrum is small (inefficient-acceleration 
case), then the energy spectral index of protons is s p ~ 2, 
and the 7r°-decay 7-ray emission shows a roughly flat vF v - 
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Fig. 4. The same as in Fig. 3, but in the wide-band en- 
ergy range fro m radio to TeV 7-rays. R adi o and X-ray data 
are ta ken from lAharonian et al.l (|2006( ) and lTakahashi et al.l 
(|2008t l. respectively. 



spectrum, vF v oc u°, in the GeV-TeV band (model Ila/IIIa 
in Table[T|). The predicted flux is marginally consistent with 
the EGRET upper limit, vF v » 5 x 10~ n erg s^cm" 2 
at 1-10 GeV (|Aharonian et alj|2006t lHartman et al.lll999f) . 
On the other hand, if the magnetic field is moderately 
strong, the synchrotron cooling effect causes steepening of 
the electron spectrum over a wide energy range — typ- 
ically s e w 3 (e.g., see § 19.3, eq. (19.16) of iLongairl 
1994). In this case, leptonic IC emission in the GeV-TeV 



band again shows a nearly flat ^F^-spectrum (model lib in 
Table [I]). Therefore, these 7-ray emission models (Ila/IIIa 
a nd lib) cannot be di stinguished. This has been discussed 
in lEllison etHI (|2007l) . where B « 60 fiG. 

In summary, it may be difficult to differentiate between 
hadronic and leptonic emission by the spectral shape of 
the GeV-to-TeV 7-ray emission of Galactic young SNRs 
like RXJ1713 (Table [T|). As shown in this paper, when 
the GeV 7-ray flux is relatively low (e.g., vF v oc v ' 5 ), 
both an efficient acceleration model with hadronic 7-ray 
emission (model la) and a leptonic, weak magnetic-field 
model with inefficient acceleration (model Illb) may give 
similar spectra l shapes. On the o ther hand, as already 
pointed out in lEllison et al.l ((2007), when the GeV emis- 
sion is relatively bright (e.g., vF v oc v ), one may not be 
able to distinguish the hadronic model in the inefficient 
case (models Ila/IIIa) from the leptonic one with a mod- 
erately strong magnetic field (model lib). This conclusion 
may, at least qualitatively, be applicable to other young 
SNRs emitting TeV g a mma-rays, such a s RX J 0852. 0—4622 
(jKatagiri et al] 120051 : lAharonian et all l2005h . Fermi will 
likely provide us with rich information on the emission 
mechanism of RXJ1713 and other young SNRs. However, 
one should only draw conclusions with great care, even in 
the Fermi era. Probably, neutri no observation with km 3 - 
class detec tors such as IceCube (Achterberg et al]|2007t ) or 
KM3NeT (jKappes fe Consortiumll2007ft will finally resolve 
the p r oblem dCrocker et al jl2002t I Alvarez-Muniz fc Halzenl 
2002t iKistler fe Beacoml 120061: lyissani fe Villantd l2008t 



I 
> 




12 13 
log (Energy / eV) 

Fig. 5. ^F„-spectra of TeV 7-rays (solid line) and fi- 
neutrinos (dashed line) calculated within the hadronic two- 
zone model in the efficient acceleration case. The solid line 
is the same as in Fig. 2. We have averaged the vacuum os- 
cillation effects among neutrinos. The dotted line shows 
the daily ave raged atmospheric neutrino flux expected 
in K M3NeT (jKappes fe Consortium! 120071 : iKappes et all 
12007ft . 



Fig- El if the observed TeV 7-ray emission is hadronic, then 
the expected neutrino spectrum at the source is above the 
atmospheric neutrino background at around 5-10 TeV, 
which may become 'the smoking gun' of proton accelera- 
tion in Galactic young SNRs. 
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